Recent Progress of the Low-Dimensional Spin- Wave Theory 
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A modified spin-wave theory is developed and applied to low-dimensional quantum magnets. 
Double-peaked specific heat for one-dimensional ferrimagnets, nuclear spin-lattice relaxation in fer- 
rimagnetic chains and clusters, and thermal behavior of Haldane-gap antiferromagnets are described 
within the scheme. Mentioning other bosonic and fermionic representations as well, we demonstrate 
that spin waves are still effective in low dimensions. 



PACS numbers: 75.10.Jm, 75.30.Ds, 75.50.Xx 
I. INTRODUCTION 

It may be spin waves that first come into our mind 
when we discuss low-energy magnetic excitations. The 
spin-wave concept— local excitations can be stabilized 
by recovering the original symmetry— is widely valid in 
physics. The ferromagnetic spin- wave theory was initi- 
ated by Bloch [1] and quantum-mechanically formulated 
by Holstein and Primakkof [2] . Dyson [3] considered cor- 
relations between spin waves [4-6]. Anderson [7] and 
Kubo [8] developed the scheme to antiferromagnets. Such 
theoretical arguments motivated numerous experiments 
[9-12]. Nowadays the spin- wave picture still stimulates 
many researchers to further explorations. 

On the other hand, we may have a vague but fixed 
idea that spin waves are no more valid in low dimensions 
where quantum as well as thermal fluctuations break 
down any magnetic long-range order. Low-dimensional 
materials were synthesized one after another, but the 
spin-wave theory was not able to play any effective 
role in studying their quantum behavior. While dy- 
namic properties can be potentially revealed in terms 
of spin waves [13,14], their application [15,16] was re- 
stricted to sufficiently low temperatures assuming the on- 
set of the three-dimensional long-range order. In such 
circumstances, (C6HnNH 3 )CuCl3 [17], which turned 
out a model compound for one-dimensional ferromag- 
nets, stimulated public interest in low-dimensional ther- 
modynamics. Numerically solving the thermodynamic 
Bethe-ansatz equations for the one-dimensional Heisen- 
berg model 
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Takahashi and Yamada [18] found that the spin-i fer- 
romagnetic specific heat and magnetic susceptibility are 
expanded as 
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at low temperatures, where * = —k^T/J (J < 0). 
Schlottmann [19] also obtained similar results. These 
findings not only settled the long-argued problem of the 
low-temperature behavior of one-dimensional ferromag- 
nets [20-24] but also motivated revisiting the subject by 
the spin- wave scheme. 

The conventional spin-wave theory [2] applied to one- 
dimensional ferromagnets gives no quantitative informa- 
tion on the susceptibility but reveals the correct lead- 
ing term of the specific heat at low temperatures. Spin 
waves still look valid in low dimensions. The difficulty 
of diverging magnetization originates in uncontrollable 
thermal excitations within the conventional spin-wave 
scheme. Then we may have an idea of constraining spin 
waves to keep zero magnetization [25,26] , which is reason- 
able for isotropic magnets. The thus-calculated thermal 
properties of spin-S* Heisenberg chains, 
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coincide very well with Eqs. (1.2) and (1.3) in the case 
of S = h. This scheme was further applied to two- 
dimensional square-lattice [27,28] and frustrated [29,30] 
antiferromagnets and is now generically referred to as the 
modified spin-wave theory. 

The spin-wave theory thus expanded toward lower di- 
mensions and higher temperatures. However, many sys- 
tems are still left for spin waves to explore. There are few 
attempts [31] at describing one-dimensional antiferro- 
magnets in terms of spin waves beyond the naivest appli- 
cation [7,8]. In the case of antiferromagnets, one and two 
dimensions are essentially distinguished. The Neel order 
exists in Heisenberg square-lattice antiferromagnets [32], 
whereas no long-range order in Hisenberg liner-chain an- 
tiferromagnets. Correspondingly, spin waves claim that 
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the ground-state sublattice magnetizations are conver- 
gent in two dimension but divergent in one dimension. 
Hence spin waves have been considered to be useless for 
one-dimensional quantum antifcrromagncts. It is also 
unfortunate that poor progress has been made in spin- 
wave treatments of one-dimensional quantum fcrrimag- 
nets, where interesting mixed magnetism [33] lies. 

Thus motivated, we first construct a modified spin- 
wave theory for one-dimensional fcrrimagnets with vari- 
ous lattice structures. The thermal and dynamic proper- 
ties are fully revealed. Then we proceed to ferrimagnetic 
clusters, which is, to my knowledge, the first application 
of spin waves to zero dimension. Finally we demonstrate 
a modified spin-wave description of spin-gapped antifer- 
romagnets and compare it with a Schwinger-boson mean- 
field approach. 



the quantity 5 — s suppresses the divergence in Eq. (2.2), 
the spin-wave scheme is better justified with increasing 
5/ s as well as Ss [36] . 

For ferrimagnets, because of no quantum divergence, 
we can consider modifying the conventional spin-wave 
scheme either according to the Takahashi scheme [26], 
that is, in diagonalizing the Hamitonian, or in the stage of 
constructing thermodynamics. Since the original Taka- 
hashi scheme ends in a poorer description of thermal 
properties, here we take the latter approach. Introducing 
bosonic operators for the spin deviation in each sublattice 
via the Holstein-Primakoff transformation 



S£ — y25 a} n a n a n , 5^ — 5 — ajja„ , 
4 = b l\J 2 « - bt,b n , s z n = -s + b\b n , 



(2.3) 



II. ONE-DIMENSIONAL QUANTUM 
FERRIMAGNETS 

Let us start our discussion from one-dimensional 
Hcisenberg ferrimagnets composed of alternating spins 
5 and s: 
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H = J ^ (S n ■ S n + S n ■ S n +l) 



(2.1) 
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When we define quantum spin reduction r as Ms — 
N(S - t) and M s = -N(s - t), where M s and M s are 
the ground-state sublattice magnetizations, spin waves 
show that [34] 
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This integral is divergent at k = in the case of 5 = s, 
which is the problem of infrared-diverging magnetization. 
Ferrimagnets, however, get rid of this difficulty. Equation 
(2.2) gives t ~ 0.305 at (5, s) = (1, |), for example, while 
the correct value turns out 0.207(1) [35]. The spin- wave 
scheme is trivially valid for large spins. Considering that 



we expand the Hamiltonian (2.1) with respect to 1/5 
(1/45 more strictly [37]) as 



U = -2NJSs + Hi+H a + OiS- 1 ) 



(2.4) 



where we have treated 5 and s as 0(5) = O(s). We 
consider the two-body term Ho of O(S ) to be a per- 
turbation to the one-body term Hi of 0(5 1 ). Via the 
Fourier transformation 
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with a being the lattice constant and then the Bogoliubov 
transformation 



,t - 



Pi 



= aj^cosh^ 
a/jSinhi^ 



- frfeSinhflfc , 
6fcCOsh0fc , 



(2.6) 



with tanh26*fc = 2V5scos(a/c)/(5 + s), we obtain a com- 
pact expression 

H t = E t + J2 [^{k)a\a k + ^(fc)4&] . (2-7) 
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FIG. 1. Dispersion relations of the ferromagnetic and antiferromagnetic elementary excitations (a) and temperature depen- 
dences of the specific heat (b) and the magnetic susceptibility (c) for the (S, s) = (1, |) ferrimagnetic Heisenberg chain. The 
linear (dotted lines) and interacting (solid lines) spin-wave calculations are compared with quantum Monte Carlo estimates. 



2 



where Ei and lof are the 0(S l ) quantum corrections to 
the dispersion relations. There exist two branches of spin 
waves in ferrimagnets, one of which reduces the ground- 
state magnetization and is thus of ferromagnetic aspect, 
while the other of which enhances the ground-state mag- 
netization and is thus of antiferromagnetic aspect. The 
ferromagnetic and antiferromagnetic spin waves convinc- 
ingly exhibit a quadratic dispersion and a gapped spec- 
trum, respectively, as is shown in Fig. 1(a). 

While the formulation is detailed in Ref. [38], we here 
mention the core idea of our ferrimagnetic modified spin- 
wave theory— how to control the number of bosons. The 
zero-magnetization constraint, which works well in ferro- 
magnets, is not relevant to antiferromagnetically coupled 
spins, keeping constant the difference between the num- 
bers of the sublattice bosons instead of the sum of them. 
Then we may consider that the staggered magnetization 
should be kept constant and introduce an alternative con- 
straint as 



^ (5 1 + s)nl 

<r=± k J{S - s) 2 + AS s sin 2 (ak) 



= (S + s)N, (2. 



where n k = (a[a k ) and n\ = {0[(ik)- Equation (2.8) 
claims that the thermal fluctuation should cancel the stag- 
gered magnetization. The thus-calculated thermal quan- 
tities are shown in Fig. 1. The ferromagnetic features 
of ferrimagnets at low temperatures, the T^-vanishing 
specific heat and the T _2 -diverging susceptibility, are re- 
produced well. The inset in Fig. 1(c) shows that the 
interacting modified spin-wave description is highly pre- 
cise at sufficiently low temperatures. As for the specific 
heat, the inclusion of the spin- wave interactions corrects 
the position of the Schottky peak, which is consistent 
with the observations of Fig. 1(a). 

There are many model compounds for the Hamiltonian 
(2.1). Kahn et al. [39] systematically synthesized spin- 
(S, i) quasi-one-dimensional ferrimagnets composed of 
two kinds of transition metals. More complicated align- 
ments of mixed spins [40,41] were also synthesized and 
investigated in terms of modified spin waves. 



III. DOUBLE-PEAKED SPECIFIC HEAT 

The ferromagnetic and antiferromagnetic excitations 
coexistent in one-dimensional ferrimagnets lead to the 
mixed features of the specific heat: the T 1 / 2 initial be- 
havior but the Schottky-type peak at mid temperatures. 
More interesting energy structure may be expected in fer- 
rimagnets. Weak magnetic field applied to ferrimagnets 
indeed causes a double-peaked specific heat [42]. When 
the antiferromagnetic gap lies far apart from the lower- 
lying ferromagnetic band, there is a possibility of such 
observations in isotropic systems without any field and 
anisotropy. In the case of alternating-spin chains illus- 
trated in Fig. 2(a), the linear spin- wave theory suggests 



the criterion as S > 2s. A series of bimetallic chain 
compounds [39] is interesting in this context, but even a 
combination of Mn (S — |) and Cu (s — \) turns out 
not enough for a double-peaked structure of the specific 
heat [43]. 

Metals of different kinds are not necessary to ferrimag- 
nets. Caneschi et al. [44] synthesized hybrid ferrimagnets 
comprising metals and organic radicals. There are fur- 
ther solutions to ferrimagnets. Figures 2(b) and 2(c) il- 
lustrate ferrimagnets of topological origin [45,46], which 
possibly exhibit a double-peaked specific heat [47]. Fig- 
ure 3 shows the modified spin-wave calculations for the 
spin-i trimeric chains depicted in Fig. 2(b) with particu- 
lar emphasis on clarifying how an extra peak grows. The 
present scheme has the advantage of visualizing each ex- 
citation mode making a distinct contribution to the ther- 
mal behavior. As the ratio J2/J1 moves away from unity, 
the system begins to divide into small clusters and every 
excitation mode becomes less dispersive. With a suffi- 
ciently small but finite ratio, an extra peak appears at 
low temperatures. The model compound Sr3Cri3(P04)4 
indeed exhibits a double-peaked specific heat satisfying 
the condition J2/J1 ~ 0.1 [45]. A double-peaked spe- 
cific heat is possible for the tetrameric chains depicted 
in Fig. 2(c) as well [47]. Extensive measurements on 
the model compound Cu(3-Clpy)2(N3)2 (3-Clpy = 3- 
chloropyridine) [46,48] and its analogs are encouraged. 

Considering the difficulty of calculating thermody- 
namic properties at very low temperatures by numerical 
tools such as quantum Monte Carlo and density-matrix 
renormalization-group methods, the modified spin-wave 
scheme can play an effective role in investigating low- 
dimensional magnets as well. We can clearly understand 
what type of fluctuations is relevant to the phenomena 
of our interest in terms of spin waves. Indeed any spin- 
wave description is trivially less quantitative in systems 
with strong bond alternation and/or frustration [49], but 
there may be an idea that quantum phase transitions can 
be detected through the breakdown of spin-wave ground 
states [29,30]. 



(a) 



(b) 



(c) 
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FIG. 2. Schematic representations of bimetallic chain com- 
pounds (a), trimeric intertwining double-chain compounds 
(b), and tetrameric bond-alternating chain compounds (c), 
where smaller and larger bullet symbols denote spins s — | 
and S > |, while solid and dashed segments mean antifer- 
romagnetic and ferromagnetic exchange couplings between 
them, respectively. 
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FIG. 3. The interacting spin-wave calculations of the dispersion relations (the upper three) and the specific heat (the lower 
three) of the trimeric chains illustrated in Fig. 2(b), where individual contributions of the ferromagnetic and antiferromagnctic 
spin waves are distinguishably shown by dotted and dashed lines, respectively, the sum of which is equal to the total drawn by 
solid lines. 



IV. NUCLEAR SPIN-LATTICE RELAXATION IN 
FERRIMAGNETS 

Next we consider observing ferrimagnets through 
nuclear spins. In the case of a bimetallic chain 
compound NiCu(pba)(H 2 0)3-2H 2 (pba = 1,3- 
propylenebis(oxamato)) synthesized by Kahn et al. [39], 
proton nuclei lying apart from electronic spins effec- 
tively work as probes to illuminate the correlation be- 
tween spins of different kinds peculiar to ferrimagnets. 
Considering the electronic-nuclear energy-conservation 
requirement, the Raman process should play a leading 
role in the nuclear spin-lattice relaxation. The Raman 
relaxation rate is generally given by 



J_ _ 47rfi.(g^B7N) 2 
Ti ~ Y.e-Ei/k*T 

^ l hf 



-Ei/k-eT 



x | (/| J2( A nS z n + B n s z n )\i)\ 2 S(E f -Ei- 7**) , (4.1) 



where A n and B n are the dipolar hyperfine coupling con- 
stants between the nuclear and electronic spins in the nth 
unit cell, = faH is the L armor frequency of the nuclei 
with 7n being the gyromagnctic ratio, and the summa- 
tion is taken over all the electronic eigenstates \i) 
with energy Ei. Taking account of the significant dif- 
ference between the electronic and nuclear energy scales 
(Hujn < 10~ 5 J), Eq. (4.1) is expressed in terms of modi- 
fied spin waves as [50] 
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J a ~" ^(Ssk) 2 + 2(5 - s)Sshuj N /J 
x [(Acosh 2 9 k - Bsmh 2 9 k ) 2 n-(n- + 1) 
+ (Asinh 2 6> fe - Bcosh 2 k ) 2 n+(n+ + 1)] , (4.2) 



where A and B are the Fourier transforms of the cou- 
pling constants on the assumption that their momentum 
dependence can be neglected [51]. 
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FIG. 4. Temperature depen- 

dences of the proton spin-lattice relaxation rate with varying 
field for the (S,s) = (1, §) ferrimagnetic model compound 
NiCu(pba)(H20)3-2H20. The interacting modified spin-wave 
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FIG. 5. Temperature dependences of the nuclear spin-lattice relaxation rate with varying field for the (S,s) = (1, |) ferri- 
magnetic Heisenberg chain calculated in terms of the interacting modified spin waves. 



In Fig. 4 the thus-calculated relaxation rate [52] is 
compared with experimental findings [51]. We have de- 
termined the coupling constants so as to well reproduce 
the observations under the condition that the probe pro- 
tons are located closer to Cu than Ni [39]. Through the 
relations A <~ dg 3 and B <~ dj 3 , where ds (d s ) is the 
average distance between the protons and the Ni (Cu) 
site in each unit, we obtain rough estimates ds — 3.8 [A] 
and d s ~ 2.6 [A], which are consistent with structural 
analyses [39]. However, the calculated field dependence 
is much weaker than the observations [51]. The field de- 
pendence can be better interpreted by the direct pro- 
cess [53] than by the Raman process. Since there are 
few nuclear-magnetic-resonance measurements on one- 
dimensional ferrimagnets, we hope that further explo- 
rations will be made from both chemical and physical 
points of view. 

Apart from particular materials, we discuss the fer- 
rimagnetic relaxation more generally When we com- 
pare Eq. (4.2) with the expression of the susceptibility- 
temperature product [38], 
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(,9MB ) 2 
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(4.3) 



we find that the ferromagnetic and antifcrromagnctic ex- 
citations interestingly make varying contribution to 1 jT\ 
according to the location of probe nuclei. Decreasing and 
increasing %T with increasing temperature are charac- 
teristic of ferro- and antifcrromagnets, respectively, and 
fcrrimagnetic %T exhibit a minimum. As long as we ob- 
serve ferrimagnetism through the susceptibility, it looks 
like a uniform mixture of ferro- and antiferromagnetism. 
On the other hand, nuclear spins can selectively extract 
ferro- and antiferromagnetic features from ferrimagnets 
according to their location. We show in Fig. 5 the re- 
laxation rate as a function of A/B = r. The nuclei lie 
closer to larger (smaller) spins S (s) for r > 1 (r < 1). 
Equation (4.2) claims that ferromagnetic (antiferromag- 
netic) spin waves can not mediate the nuclear spin relax- 
ation at all for r = tanh 2 ^ (r = coth 2 #fc). Considering 
the predominant contribution of k ~ to the integra- 
tion (4.2) at low temperatures and weak fields [52], the 



nuclear spins correlate only to antiferromagnetic (ferro- 
magnetic) spin waves at r ~ s/S (r ~ S/s). Figures 
5(a), 5(b), and 5(c) indeed present antiferromagnetic, fer- 
romagnetic, and mixed-magnetic features, respectively, 
where vanishing rates at low temperatures in Figs. 5 (b) 
and 5(c) are due to the Zeeman interaction and should 
therefore be distinguished from the intrinsic properties. 
Since antiferromagnetic spin waves are hardly excited at 
low temperatures, nuclear spins located so as to satisfy 
the condition r ~ (d s /ds) 3 — s/S exhibit extremely slow 
dynamics. 
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FIG. 6. Field dependences of the nuclear spin-lattice re- 
laxation rate with varying temperature for the (S,s) = (1, 1) 
fcrrimagnetic Heisenberg chain calculated in terms of the in- 
teracting modified spin waves, where J/ks = 121 K and r = 0. 



Another interest in Eq. (4.2) is its field dependence. 
While it is usual for the relaxation rate to depend on an 
applied field through the electronic Zeeman energy here 
is further field dependence in the prcfactor. This novel 
field dependence originates in the quadratic dispersion re- 
lation of ferrimagnets and may arise from any nonlinear 
dispersion at the band bottom in more general. In Fig. 
6 we show the relaxation rate as a function of an applied 
field, setting the parameters for NiCu(pba)(Il20)3-2H20 
[39]. At low temperatures and moderate fields, 1/Ti be- 
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haves like oc 1/yH, where exhibits a sharp peak at 
k = and thus the integration (4.2) can be approxi- 
mately replaced by the k — contribution. With increas- 
ing temperature, the predominant k — contribution 
smears and the k integration ends up with a logarithmic 
field dependence. With increasing field, these unique ob- 
servations are all masked behind the overwhelming Zee- 
man effect. The l/>/~H dependence of our interest should 
be distinguished from the diffusion-dominated dynamics 
[54-56], which appears at high temperatures originating 
from transverse spin fluctuations. We hope that ferri- 
magnetic dynamics in a field will be more and more stud- 
ied from the experimental point of view. 



V. NANOSCALE MOLECULAR MAGNETS 

Mesoscopic magnetism [57] is one of the hot topics in 
materials science, where we can observe a dimensional 
crossover on the way from molecular to bulk magnets 
[58]. Metal-ion magnetic clusters are thus interesting 
and among others is [Mn 12 Oi 2 (CH 3 COO)i 6 (H 2 0)4] [59] 
(hereafter abbreviated as Mni2), for which quantum tun- 
neling of the magnetization [60,61] was observed for the 
first time. Figure 7 illustrates the Mni 2 cluster, where 
four inner Mn 4+ spins and eight outer Mn 3+ spins are 
directed antiparallel to each other and exhibit a novel 
ground state of total spin S = 10 [62]. Assuming the sta- 
bility of the collective spin of S = 10, the Mni 2 cluster is 
often treated as a single spin-10 object [63]. However, re- 
cent electron-paramagnetic-resonance measurements [64] 
suggest a possible breakdown of the spin-10 description 
even at low temperatures. We have to consider the intr- 
acluster magnetic structure for further understanding of 
nanoscale magnets. Since NMR measurements on the 
Mn i2 cluster have recently made remarkable progress 
[65-69], we consider a microscopic interpretation [70] of 
them. 

The total spin states in the Mni 2 cluster is too large 
even for modern computers to directly handle. Indeed 
several authors [71,72] have recently succeeded in obtain- 
ing the low-lying cigenstates of the microscopic Hamil- 
tonian, yet very little is known about the intracluster 
magnetic structure. There are four types of exchange 
interactions between the Mn ions forming the Mn 12 clus- 
ter, as is illustrated in Fig. 7, but their values are still 
the subject of controversy. Taking account of the current 
argument, we consider three possible sets of parameters. 
The Hamiltonian is defined as [73] 
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"H = - ^ [2Jis„ • S n + 2J 2 (s„ ■ S„ + S n ■ s„+i) 
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+2J 3 (s n ■ s n+1 + -s n ■ s n+2 ) 
+ 2J4(S n ■ S n + S n ■ S n +\) 

+D 2 (Sz) 2 +D 3 (S*) 2 +gv B H( Sn + S z n + S z n ) , (5.1) 
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FIG. 7. Schematic plot of the Mni2 cluster. 

where s n , S n , and S n are the spin operators for the 
Mn(l) (spin § ee s), Mn(2) (spin 2 ee S), and Mn(3) 
(spin 2 = S) sites in the nth unit, respectively. The pa- 
rameter sets (a), (b), and (c) listed in Table I are based 
on the calculations in Refs. [74], [73], and [72], respec- 
tively. The early arguments (a) and (b) agree on J\ 
predominating over the rest, which is supported by re- 
cent calculations [70,71], whereas the latest claim (c) is 
completely against such a consensus, throwing doubt on 
previous investigations such as the eight-spin modeling of 
the Mni 2 cluster [75]. As for the anisotropy parameters, 
there is much less information. When the molecule is 
treated as a rigid spin-10 object, the macroscopic uniax- 
ial crystalline anisotropy parameter D is determined so 
as to fit the zero-field separation between the lowest two 
levels, A ~ 14 K [62,73]. Hence it is natural to choose 
the local single-ion anisotropy parameters D 2 and £> 3 , de- 
scribing the Jahn- Teller-distorted Mn 3+ ions [76] , within 
the same scheme. 

We introduce the bosonic operators within the lin- 



ear spin-wave scheme as s n 



-S + a n ,a n 1 



2sa' nl ; — 



2Sa 



n,2'i 



s 

S — a n 3 a n 3 , = \/2Sa n £, which can be justified in the 
low-temperature range where 55 Mn NMR measurements 
[67-69] are performed. While the same type constraint 
as Eq. (2.8) is given by 

2 2 "fc." 2 IVw(fc)| 2 -85 + 4s, (5.2) 



k n=l,2,3 



m=l,2,3 



providing we define the Bogoliubov transformation as 

«M = — ^n(fe) -ipi2(k)bl 2 +ipi3{k)bk,3, 
ak,2= ip2i(k)b kA +?A 22 (fc)&M --0 23 ( fc ) & L 3 : ( 5 - 3 ) 

Ofe,3= ^3l( fc ) & M +^32( k ) b k,2 -^33( fc ) & !,3' 

this is not applicable to the Mn i2 cluster as it is. In 
isotropic ferrimagnets, there exists a zero-energy excita- 
tion and therefore a certain number of bosons naturally 
survive at low temperatures. Once a gap A opens, which 
is the case of the anisotropic Mn i2 clusters, the boson 
number should exponentially decreases as cx e ~ A/kBT at 
low temperatures, but the constraint (5.2) still keeps it 
finite even at T — > 0. Hence we replace Eq. (5.2) by 
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FIG. 8. Semilog plots of the 55 Mn nuclear spin-lattice relaxation rates as functions of temperature under no field for the 
parameter sets listed in Table I. The linear modified spin- wave calculations, dashed [Mn(l)], dotted [Mn(2)], and solid [Mn(3)] 
lines, are compared with measurements [67], o [Mn(l)j, o [Mn(2)], and x [Mn(3)]. 
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FIG. 9. Semilog plots of the 55 Mn nuclear spin-lattice relaxation rates as functions of an applied field at T = 1.4 K for the 
parameter sets listed in Table I. The linear modified spin- wave calculations, dashed [Mn(l)], dotted [Mn(2)], and solid [Mn(3)] 
lines, are compared with measurements [67], o [Mn(l)j, o [Mn(2)], and x [Mn(3)]. 



E E »fc." E \i>mn(k)\ 2 = (8S + 4 S ) C - A / k * T . 
k n=l,2,3 m=l,2,3 

(5.4) 

This is quite natural modification of the theory, because 
the new constraint (5.4) remains the same as the autho- 
rized one (5.2) except for the sufficiently low-temperature 
region k^T < A. It is also convincing that Eq. (5.4) 
smoothly turns into Eq. (5.2) as A —* 0. 

The thus-calculated Raman relaxation rates for the 
55 Mn nuclei as functions of temperature and an applied 
field are compared with measurements [67] in Figs. 8 
and 9, where the coupling constants A4 are the only 
adjustable parameters in reproducing both temperature 
and field dependences and have been chosen as Table II. 
The parameter set (b) best agrees to the experimental 
findings. At the same time, we are surprised to learn 
that the latest argument (c) is also reasonable. Since 
the present calculations are not so quantitative as to dis- 
tinguish the two scenarios (b) and (c), the intracluster 
magnetic structure should be reexamined in all its as- 
pects. The assumption of the leading Raman process is 
less trivial for the Mn 4+ ions. A recent experiment [68] 
suggests that the hyperfine field of the Mn 3+ ions is in 
fact anisotropic and makes a predominant dipolar con- 



tribution, whereas that of the Mn 4+ ions is isotropic and 
results in the Fermi contact. In this context, it is in- 
teresting to compare carefully the theoretical (Af 1 ) and 
experimental (Af) findings for the coupling constants. 
Assuming the set (b), Af ~ 2.5 Af, Af ~ 1.7Af , and 
Af ~ 1.2A§ X . Somewhat larger deviation of the the- 
ory from the experiment for A\ implies that the nuclear 
spin-lattice relaxation on the Mn(l) site may not be Ra- 
man active primarily but be strongly influenced by the 
surrounding Mn 3+ ions. The present theory is distinct 
from the phcnomenological interpretation [65] assuming 
phonons to mediate the relaxation. We hope microscopic 
investigations of this kind will contribute toward the total 
understanding of mesoscopic magnetism. 



VI. SPIN-GAPPED ANTIFERROMAGNETS 

We have fully discussed the systems with macroscop- 
ically degenerate magnetic ground states, where the dif- 
ficulty of diverging number of bosons does not occur un- 
less temperature is sufficiently high. From this point of 
view, one-dimensional antifcrromagnets are harder for 
spin waves to treat, where the problem lies already in 
the ground state. However, here are not a few interest- 
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ing systems such as Haldane-gap antiferromagnets [77] 
and antiferromagnetic spin ladders [78]. The antiferro- 
magnetic modified spin-wave theory is not so successful 
as that for ferro- and ferrimagnets, but it can still play 
a helpful role in our explorations. We close our discus- 
sion by demonstrating bosonic treatments of integer-spin 
Hciscnbcrg chains [79]. 

We are expected to avoid the quantum divergence of 
sublattice magnetizations maintaining the core idea of 
constant bosons in number. We introduce the sublattice 
bosons as Eq. (2.3), where S — s, and diagonalize an 
effective Hamiltonian 



n = n 



2JA^«a„ 



(6.1) 



instead of TL 1 where the Lagrange multiplier A is deter- 
mined by the condition 



(6.2) 



Within the conventional spin-wave theory, spins on one 
sublattice point predominantly up, while those on the 
other predominantly down. The constraint (6.2) restores 
the sublattice symmetry. 

It is interesting to compare the modified spin-wave 
scheme with another bosonic approach. The Schwinger- 
boson representation of the spin algebra [80,81] is also a 
common language to interpret quantum magnetism and 
its application to one-dimensional systems [82,83] is mov- 
ing ahead recently. Let us describe each spin variable in 
terms of two kinds of bosons as 



St 



a^aii 



Si 



where we impose the constraint 

a| T a iT + a^aii = 2S , 



(6.3) 



(6.4) 



on the bosons. At the mean-field level, we diagonalize the 
Hamiltonian with the constraints imposed on the average 
and obtain the dispersion relation as 



uj k = 2 V / X 2 - n 2 cos 2 (afc) , 
where A and f2 are determined through 
2A(2n fe + 1) 



(6.5) 



Y \f\ 2 -n 2 cos 2 (afc) 
(2n/c + 1) cos 2 (afc) 



E 



(2S+1)L, 



L. 



(6.6) 



with rife = [ e " fc/fcBT 
is expressed as 



^\ 2 -n 2 cos 2 (afc) 

. The magnetic susceptibility 
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FIG. 10. Temperature dependences of the specific heat 
(a) and the magnetic susceptibility (b) for the S — 1 an- 
tiferromagnetic Heisenberg chain calculated by the linear 
modified spin waves (LMSW), the perturbational (PIMSW) 
and full-diagonalization (FDIMSW) interacting modified spin 
waves, the Schwinger bosons (SB), and a quantum Monte 
Carlo method (QMC). 



knT 



(6.7) 



k cr=± 



while the internal energy as 



E = JL[S 2 + 2fl 2 - 2X(2S + 1)] + 2^(2n fc + l)uj k , 

k 

(6.8) 

where we have corrected the mean-field artifact [80] of 
double counting the degrees of freedom for the Schwinger 
bosons. 

In Table III all the present findings for the ground- 
state energy and the lowest excitation gap are compared 
with quantum Monte Carlo calculations [84]. As for 
the modified spin-wave calculations, we can treat the 
O(S ) correlations in two different ways. One idea is 
the perturbational treatment of 7io to "Hi, which is re- 
ferred to as the perturbational interacting modified spin- 
wave scheme, while the other is the full diagonalization of 
Tii + 7^0 j which is referred to as the full-diagonalization 
interacting modified spin-wave scheme. On the whole, 
the bosonic languages well represent the ground-state 
energy but fail to describe the Haldane gap quantita- 
tively. However, they are still useful in understanding 
the nature of magnetic excitations. Although the mod- 
ified spin waves underestimate the absolute value of the 
gap, they reproduce the observed temperature depen- 
dence of the Haldane gap [85] better than the nonlinear 
a model [79,86]. Thus the spin-wave picture still looks 
valid for spin-gapped antiferromagnets. The Haldane- 
gap phase may be symbolically interpreted as a valence 
bond solid [87] , but such a picture is not strictly valid for 
the pure Heisenberg Hamiltonian. In the case of spin 1, 
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for example, even a likely ground state of bound crack- 
ions [88,89] moving in the valence-bond-solid background 
gives a variational energy [90] inferior to the interacting 
modified spin-wave estimates. It is also interesting that 
the full-diagonalization interacting modified spin waves 
and the mean-field Schwingcr bosons give the same esti- 
mate of the Haldane gap. The Schwingcr-boson disper- 
sion relation (6.5) indeed coincides exactly with that of 
the full-diagonalization interacting modified spin waves 
at T = 0. This is not so surprising because the Holstein- 
Primakoff bosons are obtained by replacing both a/j and 
aj| by (2S — aj^a/jj 1 / 2 in the transformation (6.3). 

On the other hand, the modified spin waves are defi- 
nitely distinguished from the Schwinger bosons in ther- 
mal calculation. Their findings for the specific heat 
and the magnetic susceptibility are compared in Fig. 
10. While the Schwinger-boson mean-field theory rapidly 
breaks down with increasing temperature, the modified 
spin waves still give useful information at finite temper- 
atures. The modified spin-wave calculations well repro- 
duce the overall behavior of the susceptibility converging 
into x/Lg 2 fiR = S(S + l)/3fcBT at high temperatures. 
We admit, however, that one of the fatal weak points in 
the modified spin-wave description of spin-gapped anti- 
ferromagnets is nonvanishing specific heat at high tem- 
peratures. It is not the case with ferro- and ferrimag- 
nets provided we introduce the Lagrange multiplier in 
constructing thermodynamics. The endlessly increas- 
ing energy with increasing temperature is because of the 
temperature-dependent energy spectrum, where the La- 
grange multiplier, playing the chemical potential, turns 
out a monotonically increasing function of temperature. 
Such a difficulty is common to all the bosonic approaches 
to spin-gapped antiferromagnets, where a fermionic lan- 
guage may be more effective [91-93]. 



VII. CONCLUDING REMARKS 

In the late 1980's, the low-dimensional spin-wave the- 
ory made an epochal progress, where a problem pend- 
ing in one-dimensional fcrromagnets was solved and two- 
dimensional antiferromagnets were discussed in the con- 
text of high-temperature superconductivity. More than 
a decade has passed and the spin- wave theory has grown 
more useful and extensive covering one and possibly zero 
dimensions. There are further interesting investigations 
in terms of spin waves for ladder [93-95] , frustrated [96] , 
and random-bond [97] systems. The present findings can 
be reconstructed through the Dyson-Maleev transforma- 
tion [3,98] instead of the Holstcin-Primakoff scheme. The 
modified spin-wave theory breaks the rotational symme- 
try anyway, while the Schwinger-boson representation 
[80-82] is rotationally invariant. 

Considering the rich harvest [26,38,70] in the ferro- and 
ferrimagnetic modified spin-wave theory, the antifcrro- 
magnetic modified spin-wave scheme [27,28,79] is less suc- 



cessful in one dimension in particular, where the ground- 
state properties and the magnetic susceptibility are well 
calculated, while the energy gap and the specific heat are 
poorly reproduced. Recently the Jordan- Wigner spin- 
less fermions defined along a snake-like path [91-93] gave 
a fine description of ladder antiferromagnets, a piece of 
which is listed in Table IV. We learn that the spin-gap is 
better reproduced by the spinless fermions. The bosonic 
and fermionic languages are complementary in our explo- 
rations of low-dimensional magnets. All these arguments 
stimulate further interest [f 00] in the spin-wave theory. 
The spin-wave scheme will be refined more and more in 
the 21st century. 
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TABLE I. Possible exchange interactions for the Mni2 
cluster in the unit of cm -1 . 





Ji 


■h 


J 3 


J 4 


(a) 


-150 


-60 


60 


30 


(b) 


-150 


-60 


-30 


30 


(c) 


-41 


-41 


3 


-8 



TABLE II. Estimates of the hyperfine interaction (?/iB7N^U 
in the unit of rad • Hz between the nuclear and electronic spins 
on the Mn(i) site. The linear modified spin-wave calculations 
are compared with experimental findings [67] . 





Mn(l) 


Mn(2) 


Mn(3) 


Theory (a) 


8.9 x 10 8 


1.8 x 10 9 


1.9 x 10 9 


Theory (b) 


2.5 x 10 8 


5.1 x 10 s 


5.4 x 10 8 


Theory (c) 


3.3 x 10 8 


5.1 x 10 8 


6.6 x 10 8 


Experiment 


1.0 x 10 8 


3.0 x 10 8 


4.5 x 10 s 



TABLE III. The ground-state energy per spin E s /L and 
the lowest excitation gap A for the S = 1 antiferro- 
magnetic Heisenberg chain calculated by the linear mod- 
ified spin waves (LMSW), the perturbational (PIMSW) 
and full-diagonalization (FDIMSW) interacting modified spin 
waves, the Schwinger bosons at the mean-field level (SB), and 
a quantum Monte Carlo method (QMC) [84]. 



EJL 



LMSW 

PIMSW 

FDIMSW 

SB 

QMC 



-1.361879 
-1.394853 
-1.394617 
-1.396148 
-1.40148H 



0.07200 
0.07853 
0.08507 
0.08507 
0.41048(6) 



TABLE IV. The ground-state energy per rung EJN and 
the lowest excitation gap A for the S — i antiferromag- 
netic Heisenberg two-leg ladder with intra- and interchain 
couplings of equal strength calculated by the spinless fermions 
within the Hartree (HSF) and Hartree-Fock (HFSF) approx- 
imations, the linear modified spin waves (LMSW), the per- 
turbational interacting modified spin waves (PIMSW), and a 
density-matrix renormalization-group method (DMRG) [99]. 



EJL 



HSF 

HFSF 

LMSW 

IMSW 

DMRG 



-0.838805 
-0.907183 
-1.101403 
-1.143750 
-1.156086 



0.500000 
0.450419 
0.100911 
0.112832 
0.504 
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